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In this article, we calculate the form factors and the coupling constant of the vertex D*D s (f> 
using the three-point QCD sum rules. We consider the contributions of the vacuum condensates 
up to dimension 7 in the operator product expansion(OPE). And all possible off-shell cases are 
considered, cj>, D s and D *, resulting in three different form factors. Then we fit the form factors into 
analytical functions and extrapolate them into time-like regions, which giving the coupling constant 
for the process. Our analysis indicates that the coupling constant for this vertex is Gds*Ds</> = 
4.12 ± ONOGeV^ 1 . The results of this work are very useful in the other phenomenological analysis. 

As an application, we calculate the coupling constant for the decay channel D J — > D s 7 and analyze 
the width of this decay with the assumption of the vector meson dominance of the intermediate 
(f>( 1020). Our final result about the decay width of this decay channel is F = 0.59 ± 0.15 keV. 

PACS numbers: 13.25.Ft; 14.40.Lb 


1 Introduction 

In relativistic heavy ion collisions J/if suppression has been recognized as an important tool to 


identify the possible phase transition to quark-gluon plasma[l|. The dissociation of J/if in the quark- 
gluon plasma due to color screening can lead to a reduction of its production. People usally explained 
this phenomenon as a process of the J/t/> absorption by 7r, p or <f> mesons in a meson-exchange 


model|2|]. And we can calculate the the absorption cross sections based on the interractions among 
the quarkonia and mesons, where the hadronic coupling constants are basic input parameters. A 
detailed knowledge of the hadronic coupling constants is of great importance in understanding the 
effects of heavy quarkonium absorptions in hadronic matter. Besides, the hadronic coupling constants 
about the heavy-light mesons can also help us understanding the final-state interacions in the heavy 


quarkonium 


decays Q]. 


Furthermore, some exotic mesons have been detected in recent years Q, which 


are beyond the usual quark-model description as qq pairs. And people interpret them as quark-gluon 
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hybrids ( qqg ), tetraquark states ( qqqq ), molecular states of two ordinary mesons, glueballs, states 
with exotic quantum numbers and many others jfj]. The form factors and coupling constants play an 
important role in understanding the nature of these exotic mesons. 


However, the strong coupling constant used in the above questions can not be explained by per¬ 
turbative theories, because the associate interactions lie in the low energy region. It is fortunate that 
the QCDSR approach can help us to solve the difficulty. The QCDSR is one of the most powerful 
non-perturbative methods, which is also independent of model parameters. In recent years, numerous 
research articles have been reported about the precise determination of the strong form factors and 


coupling constants via QCDSR, light-cone QCDSR or lattice calculation[6|-|8|]. And many strong cou¬ 
pling constants have been determined by different groups, for example, D*D S K , Ds*DI \, B*B C T, 


B*B c *Ij, B*BK , J/ipD*D s , J/i/;D s D s , J/ipD*D*, D*D s r]'^M\M- In this work, we use the QCDSR 
formalism to obtain the coupling constant of the meson vertice D*D s (f>, where the contributions of 
the vacuum condensates up to dimension 7 in the OPE are considered. The results of this work are 
very useful in these phenomenological analysis mentioned above. 


It is indicated by the BaBar collaboration that T(D*) < 1.9 MeV and 1 (D r “ D ‘ ! - ss 0.94 
However, the exact value of the decay width have yet not been determined. A more exact result can 
help us understanding the nature of the meson and testing the validity of the theoretical model. As 
an application, we also give an analysis about the decay D* —> in the end of this paper, where 

the electromagnetic coupling constant Gn*D a -y will be used. This coupling constant can be easily 
obtained, when we set Q 2 = 0 in the analytical function of coupling constant Gd*d 3 cI>(Q 2 ) in Sec.III. 



The outline of this paper is as follows. In Sec.II, we study the D*D s <j) vertices using the three-point 
QCDSR. In order to reduce the uncertainties of the result, we calculate the three-point correlation 
functions: one with the vector meson </> off-shell, another with the pseudoscalar meson D s off-shell, 
and a third one with the vector meson D* off-shell. Besides of the perturbative contribution, we also 
consider the contribution of ( qq), ( qgcr.Gq ), ( g 2 G 2 ), (/ 3 G 3 ), ( qq ) 2 and ( qq)(GG) at OPE side. In Sec. 
Ill, we present the numerical results and discussions, and Sec IV is reserved for our conclusions. 


2 QCD sum rules for the D*D s (j) 

In this work, the D*D s <f> is a vector-pseudoscalar-vector (VW) vertex. With each meson off-shell, 
we write down the three-point correlation functions: 


H tM)=i 2 f d A xd 4 ‘ye ip ' x+i ( p - p ' > v (0|T{J 5 (a:)j jU (y) J^(0)}0|) 
H %(p,p') =i 2 j d*xd*yj*'■**<*-*)■* (0\T{Mx)J 5 (y)jl (0)}0|> 

n °hp,p') = 


(1) 

( 2 ) 

( 3 ) 
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where T is the time ordered product and Jl(x), Js(x) and (:c) are the interpolating currents of the 
mesons D*, D s and <j> respectively: 


Jt(x) = s(x)j v c(c r) 

(4) 

J 5 (x) = c(x)ij 5 s(x) 

(5) 

jp.(x) = s(x) r y l j,s(x) 

(6) 


According to the QCDSR, these correlation functions can be calculated in two different ways: using 
hadron degrees of freedom, called the phenomenological side, or using quark degrees of freedom, called 
the OPE side. In the following we will obtain the sum rule according to above formulations. 

2.1 The phenomenological side 


We insert a complete set of intermediate hadronic states with the same quantum numbers as the 
current operators Jl(x), J$(x) and j^x) into the correlation functions H'^ l/ (p,p'), H^‘(p,p’) and 
(p,p') to obtain the phenomenological representations. After isolating the ground-state contribu¬ 
tions, we get the following functions for the mesons <j>, D s and D* off-shell cases. 


ytphen{4>) 


Y\phen(D s ) 


n 


phen(D *) 

pi* 


0 P 2 + ™? D ,) (q 2 + m|) {p' 2 + m 2 Da ) + 

= - - ——-l- h r 

0P 2 + ml) (q 2 + m 2 Ds ) ( p' 2 + m 2 D ,) 

-CG^(q 2 )p^e^ 

(p 2 + m 2 De ){q 2 + m 2 D ,)(p' 2 + ro|) + 


(7) 

( 8 ) 

(9) 


where C = ^ Ds mp , 3 ^ P f mP ' a and h.r. stand for the contributions of higher resonances and contin- 

(m s +m c ) ° 

uum states of each meaon. And in the derivation, we have used the following effective Lagrangian £ 
and definitions for the decay constants fD a and /^: 


£ = G DtD ^e aft xr(.d a D*+PD-d x + d a D* s ~P Dfd x )<fr T (10) 

<0|J„(0)|£>;(p)) = (11) 

(0| J 5 (0)\D s (p')) = f Ds m 2 Ds /(m s + m c ) (12) 

(0|Jm(0)I ^( q )) = ( 13 ) 


where Q, and are the polarization vectors. From Eqs.(7)~(9), we can see that there is only one 
tensor structure to work within the formalism of the QCDSR. 


2.2 The OPE side 
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Now, we briefly outline the operator product expansion for the correlation functions Eqs.(l)~(3) 
Firstly, we contract the quark fields in the correlation functions with Wich’s theorem. 

nj$ = -i 2 J d A xd A ye ip ' x+l(p - p ' )v tr{i l5 s mn {x - y)j„s nk (y - 0)j v c km {0 - x)} (14) 

IlJ™ = - i 2 J d i xd A ye ip ' x+i ^- p ' )y tr{'y v c mn {x - y)i-/ 5 s nk {y - 0) 7M s fcm (0 - re)} (15) 

n^ :) = -i 2 J d 4 xd 4 ye ip ' x+i ^- p ' )y tr{^s mn (x - y) lv c nk {y - 0)ry 5 s fem (0 - a:)} (16) 

Then, we replace the c and s quark propagators c l °{x) and s lJ (x) with the corresponding full 
propagators |l2(, 


Sij(x) 


idif! dijUi s Sij ( ss ) td'/j:/:m s ( ss ) 

2n 2 x 4 47t 2 x 4 12 48 

WsGaptijifo 00 + cr " /3 ^) i5ijX 2 tfg 2 (ss) 2 

32tt 2 :e 2 7776 

{Sj^Si) 7/i 

4 + '"’ 


t^ar ( sg s aGs) iSijX 2 ^m s ( sg s aGs ) 

192 + 1152 

<%z 4 (ss) (g 2 GG) (sja^Si) o> 

27648 8 

(17) 


Cy- (x) — 


+ 


d A ke~ lk - 


9 sG™pt™j a 010 (ft + m c ) + {ft + m c )cr a0 


(27r) 4 J ^ k — m c 4 (fc 2 — m 2 ) 2 

g s g a G^-(/ A/?a + / Aa/? ) 9 2 (t a t%G a a/ ,G b ^(f a ^ + f a ^ v + .T^) 


3 (k 2 — m 2 ) 


2^4 


4(fc 2 — m 2 ) 5 


,( 18 ) 


/ Aa/3 _ ^ + mc ) 7 A(^ + 77i c ) 7 “(|+ m^Qi + m c ) (19) 

/ a/3/iI/ = (# + m c ) 7 “(|if + m c )') 0 {ft + m c ) 7 p ((iS + m c )Y{ft + m c ) (20) 


where (g 2 GG) = (g 2 G'! l pG na0 ), t n = the A n is the Gell-Mann matrix, D a = d a — ig s G™t n , and the 
i,j are color indices. Then we compute the integrals both in the coordinate and momentum spaces, 
and obtain the correlation functions. Finally, the correlation functions can be divided into two parts: 


j-j-OPS(M) jjpert(M) _|_ jjnon— pert(M) (21) 

where M is the off-shell meson(M = (p,D S: D*). Using dispersion relations, the perturbative term for 
a given meson M off-shell can be written in the following form: 

-I poo poo pert(M) , 2\ 

np t(M W) = -^/ / p r 2 ,( , ,q l2 \ dsdu (22) 

M 4tt- J 0 J q (s -p 2 ){u -p u ) 

and the quantities s = p 2 ,u = p' 2 and q = p—p'. We put all quark lines on mass shell using Cutkosky’s 
rules(Fig.l (a) and (6)) and obtain the spectral density pft i e J t< ' M \s,u,q 2 ) 




















FIG. 1: The perturbative contributions for <j>, D s and D* off-shell. The Dashed lines denote the 

Cutkosky cuts. 


p^ D ‘\s,u,q 2 ) 


3 q 2 (s + u — q 2 + 2m 2 — 2m 2 ) Q l0 

- 7 = (m c — m s ) - -j -^- m s e^app P (23 

VA L A(s,u,g 2 ) J 

3 [ Au-q 2 )(s + u-q 2 )-2s(u + m 2 c -m 2 ) 1 , Q 

- 7 = (to c — m s )~ -—- -j-A ——A-£-— - m c e^cpP p p 

VA L A (s,u,q-) J 


f%; t{D: \s,u,q 2 ) 


3 I". . u(s + u - q 2 ) - 2u{m 2 - m 2 + u - q 2 ) 1 , 8 

— 7 = {m c — m s )— - £ - —~m s e^pp p^ 

VA L A (s,u,q 2 ) J 


where A(s, u, q 2 = (s + u — q 2 )) 2 — 4 su. As to the non-perturbative contributions, we take into account 
the contribution of (ss), ( sga.Gs }, ( g 2 G 2 ), (/ 3 G 3 ), (ss) 2 and (ss)(GG), which are showed explicitly 
in Figs 2 and 3. It should be noticed that as the consequence of the use of the double Borel transform, 
the (f> off-shell case has only the contributions of ( g 2 G 2 ) and (/ 3 G 3 )(Fig.2). Full expressions for these 
contributions of Figs 2 and 3 for D s and D* off-shell cases can be found in Appendix A,B and C, 
where the following representations will be used: 


(-l) a + b+c 7 r 2 i 

r(o)r(6)r(c)(M 1 2 ) 6 (M 2 2 ) c (M 2 ) a 


dr(r + 1) 


a+b+c— 4^.1— b—c 


( i Q 2 

(• t + 1 )m 2 

(r + l)rn 2 

(r + 1 )m§ | 

\ T M 2 + M 2 

M 2 

tM 2 

tM$ J 


(_1 )a+b+c 2 a poo 

rabc = _i_'_1_ / dr(r + l) a + b+c ~ 5 T 

mim 2 m 3 T (a)T (b)T(c)(Mf) b (M 2 ) c+1 (M 2 ) a ~ 3 J 0 { 1 


f i Q 2 

(t + l)mf 

(r + l)m| 

(r + l)m§ } 

\ T M 2 + M 2 

M 2 

tM - 2 

tM 2 J 


exp 


(27) 
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jabc 

i mim2rri3 


(-1) 


a+b+c TT 2 i 


r(a)r(6)r(c)(Af 1 2 ) f '+ 1 (M 2 2 ) c (M 2 ) a - 3 


dr(r + 1) 


a+b+c—5^_1 —b—c 


exp < — 


Q 2 


(r + 1 )m\ (t + l)r 


tM( 


M 2 


M 2 


rMj 2 


(r + 1)77 
rM , 2 


(28) 



FIG. 2: Contributions of the condensate parts ( g 2 G 2 ) and (/ 3 G 3 ) for ^ off-shell case 

2.3 The Coupling Constant and the Meson decay 

We make the change of variables p 2 —> —P 2 ,p' 2 —> —P' 2 and q 2 —» — Q 2 and perform a double Borel 
transfornrfisl to the physical as well as the OPE sides, which involves the transformation: P 2 —> M 2 
and P' 2 —> Ml, where Mi and M 2 are the Borel parameters. Then, we equate the phenomenological 
and OPE sides, invoking the quark-hadron duality from which the sum rule is obtained. 

In order to eliminate the h.r. terms from the phenomenological side in Eqs.(7)~ (9), two continuum 
threshold parameters so and u 0 in the OPE side are introduced. These parameters fulfill the following 
relations: m 2 < So < m' 2 and ?n 2 < uq < m' 2 , where rrii and m 0 are the masses of the incoming and 
out-coming mesons respectively and m! is the mass of the first excited state of these mesons. After 
these performaions, the form factors can be written as: 
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(a) (b) (c) (d) (e) 




(u) 


(v) 


(w) 


(X) 


(y) 



FIG. 3: Contributions of the non-perturbative parts for D S (D*) off-slrell case 


GtnJQ 2 ) = 


~ 4-7T 2 Af 2 Af 2 C C (*> ^dsdu + SSalIinon-pertW] 


W+^)MfMi e 


./Ml 


JMl 


> (29) 



G 


D s 

D*D S 


= 


~ 4n 2 M 2 M 2 


11° C P pert(Ds) ^ u ’Q 2 ) e *** e “ 2 dsdu + 


non—pert(D s 


C' -ml/M 2 -ml /Mf 

( Q 2 + m 2 Ds )M 2 M 2t e 


,(30) 


Gd\d s AQ 2 ) — 


-4WFCCr t(D:) (s,«,Q 2 )e ^ dsdu + ^[IT“»-P ert ^)] 


C , Q— m n- /Ml - m 2 /M 2 

(Q 2 +m 2 D ,)M 2 M 2 ^ e 


(31) 


where ] stands for the double Borel transform. Now, we can calculate the form factors in 

the space-like region according to these above equations. However, in order to obtain the coupling 
constants, it is necessary to extrapolate these results into physical regions(Q 2 < 0 ), which is realized 
by fit the form factors into suitable analytical functions. It is indicated that we should get the 
same values for the coupling constants G^i D r4 and G^^^Q, when we take Q 2 = — m 2 , 

Q 2 = — m 2 Da and Q 2 = —m 2 D * separately. This above procedure is used to minimize the uncertainties 
in the calculation of the coupling constant, which will be quite clear in the following section. 


With the assumption of the vector meson dominance^ 1020)), the radiative decays D* —> D s 7 can 
be described by the following electromagnetic lagrangian £', 


£’ = -eQ s s 7 M sH M 


(32) 


where the A^, Q s are the electromagnetic field and the charge number. From the lagrangian £', we 
can obtain the decay amplitude [15]. 


(D s (ph(q,e)\D*(p',£)} 

= (7(9, e)\<!>{q, v)) — 2 ( D s(p)Hq,v)\D* s {p',0) 

= {Ds(p)(j>(q,v)\D*(p',^)) —g -f^m 4 ,eQ s (~i)£*i/ 

q -m^ M 

= G D * Ds - t e a P XT p' a q\!;i3r]*—^—^—? Um^eQ s (-i)e*rf (33) 

The parameters Gd*b s7 and are the coupling constant and the weak decay constant, respectively. 
p' a and q\ are the four momenta of the D s and 7 . rj £* and £7 are the polarization vectors of the <f>, 
7 and D*, respectively. The strong coupling constant Gn^Ds-y can be related to the effective coupling 
constant in the heavy quark effective Lagrangian by Eq.(10) in this paper. 

3 Results and Discussions 


Present section is devoted to the numerical analysis of the sum rules for the coupling constants. The 
decay constants and hadronic parameters used in this work are taken as = 0.229±0.003 16j, /d = 


0.257±0.006[3, f D . = 0.301±0.013[l6j, m$ = 1.019±0.020GeF[l6|, m Ds = 1.968±0.00032GeF|l6|, 
m D , = 2.112 ±0.0005GeI/Q] . The vacuum condensates are taken to be the standard values < ss >= 
— (0.8 ± 0.1) x (0.24 ± 0.01GeF) 3 Q, Q, < sg s aGs >= m 2 0 < ss >Q, Q, «o = (°- 8 ± 0.1 )GeV 2 , 
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-EL 


< g 2 s GG >= (0.022 ± 0.004)GeU 4 [18||, < f 3 G 3 >= (8.8 ± 5.5 )GeV 2 < g 2 GG And we also 

take the masses of quark m c = (1.275 ± 0.025)GeU, m a = 0.095 ± 0.005GeU from the Particle Data 
Group Q. The continuum parameters, Sq and uq in Eqs.(29)~(31), are defined as Sq = (to, + A i) 2 
and no = ( m 0 + A 0 ) , where the quantities A and A 0 are determined imposing the most stable 
Borel window. In order to include the pole and to exclude the h.r. contributions for the cases of </>, 
D s and D* mesons off-shell, the values for A^, Ad s and A£>. can not be far from the experimental 
value of the distance between the pole and the first excited state[13]. In addtion, the results of 


the form factors in Eqs.(29)~(31) should also not depend on the Borel parameters M 2 and M 2 . 
Therefore, we have to work in a region where the approximations made are supposedly acceptable 


and where the results depend only moderately on the Borel variables [pi. Using the Borel region 
5.0 < M 2 < 7.0 GeV 2 and 5.0 < Mf < 7.0GeU 2 (Q 2 = 3.0GeU 2 for off-shell), 6.0 < M 2 < 8.0 GeV 2 
and 6.0 < M| < 8.0GeU 2 ( Q 2 = l.OGeU 2 for D s and D* off-shell) we found a good stability with 
A 0 = A Ds = A D . = 0.5GeU(Fig.4). 


From the figure, we can see that the values are rather stalbe with variations of the Borel parameters, 
it is reliable to extract the form factors. Besides of the pertubative term, we can also see that (ss) 
give a considerable contribution for D s and D* off-slrell cases(Fig.4 (c)~(f)). And the contributions 
of the other condensate terms are small(< 1%). To the case of <j> off-shell, condensate parts ( g 2 G 2 ) 
and (/ 3 G 3 ) make up 1% ~ 2% of the total contributions. It should be noticed that although these 
condensates terms, all except for the perturbative term and (ss), give small contributions to the form 
factors, they have a significant influence on the following analytical functions(Eqs.(34) ~ (36)), which 
are obtained by numerical fitting. Thus, these condensates contributions should not be neglected in 
the calculation. 

The form factors G^l and G^; D ^ are shown explicitly in Fig.5 and are fitted into the 

folowing analytical functions: 


D* a D a 

D B 

D* a D s 

D* a 

D*D S 


= Aexp(—BQ 2 ) 

C 


1 + DQ 2 

C 

1 + D'Q 2 


exp(—EQ), 


exp{—E'Q ) 


where 


A = 2.964 ± 0.089GeU _1 , B = 0.1621 ± 0.0077GeV~ 2 , 

C = 2.755 ± 0.008, D = -0.1944 ± 0.0186, E = 0.256 ± 0.0265, 

C' = 2.825 ± 0.012GeU _1 , d = -0.1855 ± 0.0171GeU- 2 , E' = 0.2593 ± 0.0257, 


(34) 

(35) 

(36) 


Considering uncertainties of all the input parameters, such as quark and mesons masses, decay 
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(a) 


(b) 




\ 

0 

g. 

t 


o 




(e) 


(0 


FIG. 4: The contributions of different condensate terms in the OPE with variations of the Borel 
parameters and M| for 0((a),(b)), _D s ((c),(d)) and Z?*((e),(f)) off-shell, where A-H denote the 
perturbative term, (g 2 G 2 ), ( f 3 G 3 ), (ss), ( sga.Gs ), (ss) 2 , ( ss)(GG ) and Total contributions. 
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(a) (b) (c) 

FIG. 5: Form factors of the D*D s (f> vertex for cf> off-shell(a), D s off-shell(b) and D* off-shell(c). 

constants and the vlues of different condensates terms, we plot the upper and lower bounds of the 
form factors in Fig.5. We can see that although the uncertainties are large(about 20% ~ 30% of the 
central values), the fitted functions can reproduce the central values of the form factors well. Thus, it 
is reliable for us to extrapolate the Q 2 to the physical region Q 2 < 0 for (j>, D s and D* off-shell cases 
to obtain the coupling constant for the vertex D*D s cj). Substituting Q 2 = — m^, Q 2 = —and 
Q 2 = —m 2 D , separately in Eqs.(34) ~ (36), we obtain the values for G^ D , D ^, Gg» Ds( £ and 


GL^^S.SliO.llGeV- 1 (37) 

Gg^ = 4.24±0.47Gey- 1 (38) 

G^Zds* = 4-61 ± 0.51GeE _1 (39) 

Although the values for each off-shell case are different, they are roughly compatible when the un¬ 
certainties are taken into account, where the uncertainties in Eqs.(37) ~ (39) originate from the 
uncertainties of the fitted parameters SA , SB, SC, 5D, SE, SC', SD' and 8E' . Taking the mean values 
between the numbers presented above, we obtain the strong coupling constant for Gds*Ds 4 >'- 

G DslfDs ^ = 4.12 ± OJOGeV- 1 (40) 

In Reference [l9| . Z.G.Wang studied the D*DV vertex with the light-cone QCD sum rules. Most of 
the results in this work is analyzed to be much smaller than others (lfl| . And the coupling constant 
Gds*Ds4> is estimated to be about 0.82 ± 0.16GeE -1 which is also much smaller than our result. This 
difference is most probably due to the different input parameters and the different methods employed. 
In Reference [loj], it is indicated that the value of Gds*Ds4> is 4.07 ± 0.71GeE -1 in the framework of 
the three-point QCD sum rules. Besides of the perturbative part, the contributions of quark-quark, 
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gluon-gluon, and quark-gluon condensate are considered in this work. It is clearly that our result is 
compatible well with that of Referencjioj], which indicates to some extent the reliability of our result. 

As to the coupling constant Gd*d s 7 in the decay D* —> D s "f in Eq.(2), we can easily obtain its 
value by setting Q 2 = 0 in the analytical function(Eq.(34)): 

G flstDs7 = 2.96 ± 0.09GeR _1 (41) 


Now, it is time for us to give an analysis of electromagnetic decay D* —> D a 7 . As to its decay width, 
it can be written as the following representation: 


r = 


1 


2J+ 1 




(42) 


P = 


y/(M? - ( M 2 + ro 2 ) 2 )(M f - (M 2 - m 2 ) 2 ) 

2 Mi 


where the i and / denote the initial and final state mesons, respectively, the J is the total angular 
momentum of the initial meson, the ^ denotes the summation of all the polarization vectors, and the 
T denotes the scattering amplitudes. 


With the Eqs.(33) and (42), the decay width of D* —> D s "f can be expressed as: 

1 3 r M 2 ,+M 2 -" 


T — aG 2 D * Dal 


—Qs 

L TO0 


\P\- 


D s 


J 3M 2 , L 


2 M, 


di 


(43) 


with Q s = i, a = yA. Considering all the uncertainties of the input parameters, we finally get the 
decay width of the process D* —> Dg'y: 


T = 0.59 ± O.lbkeV 


It is indicated by the Babar collaboration that the decay width of r (D*) < 1.9MeV. And the ratio 
of the decay channel D* —> D s 7 is about 94.2% of the total width. This means that our result is 
compatible with the experimental data. Besides, Donald et al. predicted the value of this decay 
channel is T(D* —> D s 7 ) = 0.066 ±0.026^614 in Full Lattice QCD 201, which is much smaller than our 
result. Although these results are all compatible with the experimental data, it needs to be further 
testified by more experiments and theoretical calculations because of this difference. 


4 Conclusions 

In this article, we have calculated the form factors Gpl Da( p and in the space-like 

regions with (f>, D s and D* off-shell cases by three different QCD sum rules. Then we fit the form 
factors into analytical functions, extrapolated them into the time-like regions, and obtained the strong 
coupling constant Gds*Ds 4 >- This procedure help us to reduce the errors related to the method, leading 
to compatible coupling constants, as seen Eqs.(37) ~ (39). In addition, we also obtained the coupling 
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constant Gds*Ds~/ with the analytical function. With this coupling constant, we calculated the decay 
width of the electromagnetic decay D* —> D s 7 and compared our result with those of other groups. 
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Appendix A: Full expressions of the < g 2 G 2 > and < f 3 G 3 > contributions 
for d> off-shell case. 
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Appendix B: Full expressions about the condensate terms for D s off-shell 
case. 
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Appendix C: Full expressions about the condensate terms for D* off-shell 
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